Fermi condensates for dynamic imaging of electro-magnetic fields 
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Ultracold gases provide micrometer size atomic samples whose sensitivity to ex- 
ternal fields may be exploited in sensor applications. Bose-Einstein condensates 
of atomic gases have been demonstrated to perform excellently as magnetic field 
sensors |lj in atom chip [21 [3] experiments. As such, they offer a combination of res- 
olution and sensitivity presently unattainable by other methods [1] . Here we propose 
that condensates of Fermionic atoms can be used for non-invasive sensing of time- 
dependent and static magnetic and electric fields, by utilizing the tunable energy gap 
in the excitation spectrum as a frequency filter. Perturbations of the gas by the field 
create both collective excitations and quasiparticles. Excitation of quasiparticles re- 
quires the frequency of the perturbation to exceed the energy gap. Thus, by tuning 
the gap, the frequencies of the field may be selectively monitored from the amount of 
quasiparticles which is measurable for instance by RF-spectroscopy. We analyse the 
proposed method by calculating the density-density susceptibility, i.e. the dynamic 
structure factor, of the gas. We discuss the sensitivity and spatial resolution of the 
method which may, with advanced techniques for quasiparticle observation [5] , be in 
the half a micron scale. 

The density of an ultracold alkali gas is sensitive to spatially varying magnetic fields due 
to the Zeeman effect. This is the principle behind magnetic trapping: atoms in low field 
seeking states are trapped at the minima of the field. On the other hand, it can be used for 
sensing since magnetic perturbations leave marks on the density of the gas. Such magnetic 
field imaging has been experimentally demonstrated with Bose-Einstein condensates in mi- 
crotraps [H S]. The microkelvin sample of atoms is magnetically trapped at about 5-200 pm 
distance from the room temperature chip surface. Any additional perturbing magnetic field 
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5B will displace the center of the trapping potential and this can be measured by absorption 
imaging, either in situ or after ballistic expansion. The change in the trapping potential V 
is directly proportional to the additional field, i.e. 5V oc 5B. Similarly, electric fields can 
be sensed using the Stark effect, 5V oc 5E [4j. The principle of using density perturba- 
tions of an ultracold atomic gas for sensing can be extended, for instance, to gases that are 
trapped optically, that are not on atom chips but brought in the vicinity of the sample by 
other means, or that consist of fermionic atoms instead of bosons. We propose to utilize the 
pairing gap present in a fermionic superfluid for temporally and spatially resolved imaging 
of magnetic (or electric) fields. Superfluids of Fermi gases have been recently observed, for 
reviews see e.g. [HI [7j. Also degenerate Fermi gases in microtraps have been realized [8]. 

In the proposed method, the Fermi condensate is trapped, magnetically or optically, near 
the sample of interest. Magnetic fields in the sample, generated for instance by electric 
currents or even spin, cause density perturbations to the Fermi gas. The perturbations, 
providing energy and momentum to the gas, lead to collective or quasiparticle excitations. 
The sensing is initiated by having a high value for the excitation gap A. Only frequencies 
above 2A will be able to break pairs. The gap A can be controlled with a Feshbach res- 
onance (or by changing the density). Gradually changing the gap allows the isolation of 
individual frequencies: every time 2A crosses a frequency present in the magnetic field, the 
measured amount of quasiparticle excitations increases abruptly, see Figures [T] and [3] The 
quasiparticles can be detected by RF-spectroscopy [U [TUl HU IT^ . 

For spatial imaging of static fields, the following variant of the method can be used: The 
spatial dependence of the static perturbation provides momenta for the gas but no energy. 
Energy is given by modulating the gas uniformly in space, with a frequency corresponding 
to the pair breaking. In other words, the static perturbations serve as nucleation centers for 
quasiparticles under time-periodic modulation. 

Within linear response, the density response is 



where we calculate the susceptibility x with the generalized random phase approximation, 
following [13]. We solve x(q, tu) numerically from the most general form given in [13], 
without making the approximation of weak coupling strength. For the equation of state, see 
Methods. The magnetic field is taken to be of the form B = Ai6{u — Ui)ipi{q), where Ai 
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is the amplitude. The momentum part, fi{q), is due to the geometry of the perturbation 
and we assume it is independent of frequency. Then 

Sp{q, uj) = x(q, ^)<^(q) AS{uJ - uJi). (2) 

i 

All the relevant information is embedded in or rather its imaginary part, the dy- 

namic structure factor: ^(q, cu) = — l/7rlmx(q, u;). The dynamic structure factor has two 
parts: Anderson-Bogoliubov (AB) phonon which is a collective mode with frequency below 
2A, and quasiparticle excitations with frequencies above 2A, see Figure |2] The results are 
in qualitative agreement with those in [131 CH 113 IIHI [HI HI [IS]- 

The strong dependence of the qualitative behaviour of the dynamic structure factor on 
momentum. Figure [2| allows to focus on perturbations of a chosen length scale. The AB- 
phonon, or the collective modes of a harmonically trapped gas, may be used for detecting 
spatially large scale perturbations. Here we concentrate on perturbations of small size (1/g ~ 
l/(2kF)) which cause a strong quasiparticle response near and above the pair breaking 
frequencies. For sizes smaller than l/(2/cp') the quasiparticle threshold loses its dependence 
on A and approaches the free particle dispersion uj (x q^. 

In Figure [3] we show the response for different values of the gap A, in a case where the 
perturbation contains four different frequencies, with Ai = 1 for all (see Eq. ([2])). The 
response is the sum of dynamic structure factors for the four frequencies. The frequencies 
show up as prominent features in the amount of quasiparticles when the gap is varied. Note 
that for both momenta {QAkp, Ikp), the peaks caused by quasiparticle formation are very 
similar. Thus for a realistic perturbation geometry, whose Fourier transform contains several 
momenta, the signal should still be well resolved as long as the perturbation is roughly of 
the size 1/kp. 

The amount of quasiparticles can be measured by applying RF-pulse(s) at zero and/or 
negative detunings (see Methods). The gases are typically confined by a harmonic potential, 
therefore the density and the gap are not uniform throughout the gas. Figure |4] shows that 
the threshold type behaviour disappears when the trapping potential has been taken into 
account by local density approximation (see Methods), but the frequency of the perturbation 
is still visible as a maximum. However, we found that such smoothened response allows to 
isolate only a few, not very closely spaced, frequencies, unlike in the homogenous case. 

With tomographic techniques [5l [20], the RF-spectroscopy can be spatially resolved in 
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Figure 1: A Fermi condensate is trapped near the sample of interest. A gap opens around the 
Fermi level of the superfluid and sets the minimum energy of single particle excitations to the 
value of the order parameter, A. Magnetic fields with certain frequency and location in the sample 
cause density perturbations in the condensate. Only if the frequency exceeds 2A, quasiparticles 
are created, which allows sensing perturbations of different frequencies by tuning the gap. 
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Figure 2: Dynamic structure factor S* as a function of frequency, with two momenta, q = O.lkp 
(solid line) and q = OAkp (dashed line). The smaller momentum case shows the AB phonon as 
a clear peak, and the quasiparticle continuum above uj = 2A ~ 0.13, whereas when q approaches 
kp the phonon merges with the quasiparticle continuum. Both curves are scaled to unity for 
readability. 
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Figure 3: Dynamic structure factor as a function of the pairing gap, summed for four frequencies, 
0.03, 0.06, 0.09, and 0.12, with two different momenta, 0.4^^ (left) and Ikp (right). This corre- 
sponds to the amount of quasiparticles caused by a perturbing field with these four frequencies. 
The solid line shows the data with the AB phonon suppressed: for the momenta considered, the 
dynamic structure factor mainly corresponds to quasiparticle creation as the solid line closely fol- 
lows the full result. The dashed line shows the result at a finite temperature, T = O.OlEp, which 
is of the order O.bTc for the A considered. 
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Figure 4: The dynamic structure factor for a single frequency co = 0.03 and q = Ikp , averaged 



for a harmonic confinement using LDA, as a function of the gap at the center of the trap Aq. 



Since the pairing gap is always small at the edges of the trap, there is a finite response already 
for Aq > ujq/2. Decreasing Aq allows quasiparticle creation in larger areas in the trap, increasing 



the response, but once the quasiparticles can be formed also at the center of the trap, such growth 



stops. This leads to a maximum of the response at Aq = uJo/2, shown by the vertical dashed line. 
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three dimensions. In our proposed method, spatially resolved RF-spectroscopy could be 
used for accurate determination (without smoothening by the trap-averaging) of the pertur- 
bation frequencies and, naturally, for resolving the perturbation spatially (also in the static 
version of the method). Furthermore, the non- uniform density profile of the trapped gas 
simultaneously provides experiments with different gap values, which could be utilized when 
the perturbation is, e.g., a long thin wire. 

The frequencies that can at the present be resolved with the method are limited by the 
experimentally demonstrated gap values to the order of 10 kHz. At the unitarity limit, 
the gap becomes proportional to the Fermi energy, thereby higher particle numbers allow 
higher frequencies. Within linear response, which is proportional to the sensitivity is 

basically given by the time available for the measurement. We estimate the sensitivity to be 
10~® — 10^^^ Tesla (see Methods). To detect a single spin, the maximum distance of the gas 
from the surface is estimated to be about ~ 0.5 pm which is not possible due to noise and 
heating of the gas for samples at room temperature [2] but may be for those at cryogenic 
temperatures f2T[ [22] or for ones utlilizing photonic band gap materials . Using Feshbach 
resonances, the gas of Fermions can be also converted into a Bose-Einstein condensate of 
molecules [211 [25]. Thereby, a setup used for the Fermi condensate sensor proposed here 
could be easily turned into one that functions as the Bose-Einstein condensate sensor [H [3] 
as well, only with double mass of the particles which increases the sensitivity. 

At the present, several other systems than BCS-type superfluids are being pursued with 
ultracold gases: the proposed method could be extended to other gapped systems and 
thereby to new frequency ranges. This could also allow higher spatial resolution: Quite 
naturally, the spatial resolution of a response that involves interparticle correlations is 
given by the interparticle distance 1/kp, which for typical trapped Fermi gases is about 
l/{2kp) ~ 0.5/im as discussed above. It can, however, be smaller in optical lattices [6l [26] 
and, especially, the self-assembled crystals of ultracold polar molecules proposed in [27] could 
offer interparticle distances and thus resolutions in the nanometer scale. 

In summary, we have proposed to use an ultracold Fermi gas in a gapped state as a sensor 
for time-dependent and static magnetic fields. The tunable gap works as a frequency filter, 
and the locations of the perturbation act as nucleation centers for quasiparticles measurable 
with RF-spectroscopy. 
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METHODS 

We assume a two-component (pseudospins t and |) Fermi gas in a superfluid state de- 
scribed by the standard Bardeen-Cooper-Schrieffer (BCS) theory, given by the Hamiltonian 

H = ^(ek - /x)(c{kCtk + clkCik) + Ac{kc{_k + Aq_kCtk- (3) 

k 

The order parameter A and the chemical potential fi are obtained by iteratively solving the 
self-consistent crossover equations 

kc 

J E, ^ ' 



and 

kc 

3 f fk'^-fi 
2 



kc 2 

/ (^Er^^""^^^'^ - 1) + 1) (5) 





where = a/ (Z;;^ — /i)^ + is the BCS quasiparticle dispersion, np{e) = 1/(1 + e^^) is the 
Fermi function, kc is the cut-off, and kpa is the dimensionless coupling constant. 

We use interactions parameters in range > kpa > —0.66, resulting in pairing gaps 
A up to O.lEp. All our calculations are at zero temperature except the dashed line in 
Figure [3j We the maximum used kpa = —0.66 which is in the BSC limit just in order to 
be able to do the finite temperature calculation within simple BCS theory. The method 
itself is by no means limited to weak interactions, and actually all the estimates about the 
performance (frequencies, sensitivities, etc.) are done assuming that the experiments are 
done at the unitarity limit. Note that while the AB phonon is a signature of superfluidity, the 
quasiparticle creation does not require a superfluid. Therefore a gas at temperatures above 
Tc but having a pseudogap [28] could serve as well but the response would be smoothened 
due to the lack of sharp features in the density of states |T71 1^ . 

To detect the quasiparticles, RF pulses transferring atoms in one of the components | or 
I to a third internal state are applied with zero and/or negative detunings (or positive if 
there are strong Hartree contributions [5]), avoiding detunings which would break pairs. In 
this way only the quasiparticles produced by the magnetic field perturbation are observed. 
Note that the RF pulse length can be rather short, increasing the operation speed of the 
sensor, since high energy resolution is not required; actually it can be an advantage if the 
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pulse samples several negative/zero detunings simultaneously via the large linewidth. The 
quasiparticle response could be calibrated by experiments with known perturbations, e.g. 
microfabricated current carrying structures. Moreover, the static structure factor of a Fermi 
gas can be measured by Bragg spectroscopy [I9] which is also be useful for calibration. 

In order to account for the effects caused by the harmonic trapping, we have used the 
local density approximation (LDA) to average the signal over the trap. One defines a local 
chemical potential 

/i(r) = fiQ- ^mt^V, (6) 

where /xq is the chemical potential at the center of the trap, and calculates S{^{r)) at distance 
r as for a uniform system. The result is given by 



5'lda oc / S{fi{r))r^ dr 





2 \ mcu'^ 



(7) 



where fi is in the units of ^q. Note that this reasoning assumes perturbations spanning 
the whole gas. One or a few localized centers would again give sharp response, without 
the need for such trap-averaging, however, there would be ambiguity in determination of 
uj if the location of the center is not resolved too. The final state momentum-resolved 
RF-spectroscopy [29] could be useful in this context. 

When estimating the time available for the experiment, one should consider not only the 
lifetime of the gas which can be easily 100 ms - 1 s or even longer, but also the diffusion 
time of quasiparticles if high spatial resolution is aimed at. According to the measurements 
in [20], no significant diffusion happened during 5 ms. Therefore we take 10 ms and 1 s as 
the lower and upper bounds for the time available when estimating the sensitivity. 

The probability for producing an excitation with potential energy V applied for duration 
r is proportional to iVr/hl"^. Assuming the probability needed for a detectable signal (min- 
imum number of excited particles) is at best 0.01 and at worst 1, the minimum potential 
energy V is between 0.1/i/r and h/r. With the 10 ms -Is time scales given above the 
potential energy sensitivity lies between h ■ 0.1 Hz and h ■ 100 Hz = hu. 

The potential experienced by a neutral atom in the hyperfine state mp is = gfiBfnpB, 
where Hq is the vacuum permeabihty, /xb = eh /2m is the Bohr magneton, and g ^ 2 is the 
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Lande factor. Therefore, assuming that the potential energy sensitivity is hu, the magnetic 
field sensitivity is hv /2^b^f ~ 3.6 ■ 10^^^ T/Hz for mp = 9/2. With the limits for u given 
above, the sensitivity is between 10~^^ and 10^® T. 

Detection of a single spin is in principle possible. The magnitude of the magnetic field 
due to the spin of an electron is approximately B{r) = fiofiBd / 4:Trr^ , where /iq is the vacuum 
permeability, and r is the distance from the electron. Therefore the required sensitivity to 
be able to detect a single spin has an upper bound of 

^y^[Wnp_ (8) 

Conversely, assuming the potential energy sensitivity of /i ■ 1 Hz (from our estimated 

range of 0.1 - 100 Hz), the maximum distance at which the detection is possible is 

1 

(/io/UB"^F/^7rHz) 3 ^ 0.6 pm for mp = 9/2. 

ACKNOWLEDGEMENTS 

We thank J. Hecker Denschlag for useful discussions. This work was supported by the 
National Graduate School in Materials Physics, Ellen and Artturi Nyyssonen foundation, 
Academy of Finland (Project Nos. 213362, 217045, 217041, 217043) and conducted as a 
part of a EURYI scheme award. See www.esf.org/euryi. 



* Electronic address: |paivi.torma@hut.fi| 
[1] S. Wildermuth, S. HofFerberth, I. Lesanovsky, E. Haller, L. Andersson, S. Groth, I. Bar-Joseph, 

P. Kr uger, and J. Schmiedmayer, Nature 435, 440 (2005). 
[2] R. Folman, P. Kriiger, J. Schmiedmayer, J. Denschlag, and C. Henkel, Adv. At. Mol. Opt. 

Phys. 48, 263 (2002). 
[3] J. Fortagh and C. Zimmermann, Rev. Mod. Phys. 79, 235 (2007). 

[4] S. Wildermuth, S. HofFerberth, I. Lesanovsky, S. Groth, P. Kriiger, and J. Schmiedmayer, 

Appl. Phys. Lett. 88, 264103 (2006). 
[5] A. Schirotzek, Y.-L Shin, C. H. Schunck, and W. Ketterle (2008), arXiv:0808.0026. 
[6] L Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008). 
[7] S. Giorgini, L. Pitaevskh, and S. Stringari (2007), arXiv:0706.3360. 



10 



[8] S. Aubin, S. Myrskog, M. Extavour, L. LeBlanc, D. McKay, A. Stummer, and J. Thywissen, 
Nature Physics 2, 384 (2006). 



[9 

[lo: 
[11 

[12: 

[13 
[14 

[15: 

[16 

[17: 
[is: 

[19 
[20 

[21 

[22: 
[23 

[24 

[25: 
[26 

[27: 

[28: 

[29 



C. A. Regal and D. S. Jin, Phys. Rev. Lett. 90, 230404 (2003). 

S. Gupta, Z. Hadzibabic, M.W.Zwicrlcin, C. Stan, K. Dieckmann, C. Schunck, E. vanKempen, 

B. Verhaar, and W. Ketterle, Science 300, 1723 (2003). 

C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, J. H. Denschlag, and R. Grimm, 
Science 305, 1128 (2004). 

J. Kinnunen, M. Rodriguez, and P. Torma, Science 305, 1131 (2004). 

R. Cote and A. Griffin, Phys. Rev. B 48, 10404 (1993). 

A. Minguzzi, G. Ferrari, and Y. Castin, Eur. Phys. J. D 17, 49 (2001). 

H. P. Biichler, P. Zoller, and W. Zwerger, Phys. Rev. Lett 93, 080401 (2004). 

R. Combescot, S. Giorgini, and S. Stringari, Europhys. Lett. 75, 695 (2006). 

G. M. Bruun and G. Baym, Phys. Rev. A 74, 033623 (2006). 

K. J. ClialUs, R. J. Ballagh, and C. W. Gardiner, Pliys. Rev. Lett. 98, 093002 (2007). 

G. Veeravalh, E. Kuhnle, P. Dyke, and C. J. Vale (2008), arXiv:0809.2145. 

Y. Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, Phys. Rev. Lett. 99, 090403 (2007). 
J. Verdu, H. Zoubi, C. Koller, J. Majer, H. Ritsch, and J. Schmiedmayer (2008), 
arXiv:0809.2552. 

V. Dikovsky, V. Sokolovsky, B. Zhang, C. Henkel, and R. Folman (2008), arXiv:0808.1897. 
J. Bravo-Abad, M. Ibanescu, J. D. Joannopoulos, and M. Soljacic, Physical Review A 74, 
053619 (2006). 

S. Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, S. Riedl, C. Chin, J. H. Denschlag, and 
R. Grimm, Science 302, 2101 (2003). 

C. Regal, M. Greiner, and D. Jin, Nature 426, 537 (2003). 

M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. S. De, and U. Sen, Adv. Phys. 56, 
243 (2007). 

H. P. Biichler, E. Demler, M. Lukin, A. MicheH, N. Prokof 'ev, G. Pupillo, and P. ZoUer, Phys. 
Rev. Lett. 98, 060404 (2007). 

Q. Chen, J. Stajic, S. Tan, and K. Levin, Phys. Rep. 412, 1 (2005). 
J. T. Stewart, J. P. Gaebler, , and D. S. Jin, Nature 454, 744 (2008). 



